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We use kinematic data from three new, nearby, extremely low-luminosity Milky Way dwarf galax- 
ies (Ursa Major II, Willman 1, and Coma Berenices) to constrain the properties of their dark matter 
halos, and from these make predictions for the 7-ray flux from annihilation of dark matter particles 
in these halos. We show that these ~ 10 3 Lq dwarfs are the most dark matter dominated galaxies 
in the Universe, with total masses within 100 pc in excess of 10 6 Mq. Coupled with their rela- 
tive proximity, their large masses imply that they should have mean 7-ray fluxes comparable to or 
greater than any other known satellite galaxy of the Milky Way. Our results are robust to both 
variations of the inner slope of the density profile and the effect of tidal interactions. The fluxes 
could be boosted by up to two orders of magnitude if we include the density enhancements caused 
by surviving dark matter substructure. 

PACS numbers: 95.35.+d,14.80.Ly,98.35.Gi,98.62.Gq 



I. INTRODUCTION 



The census of the Local Group has changed dramat- 
ically in the last few years. Prior to the turn of the 
century, there were only eleven known satellite galaxies 
of the Milky Way (MW), with a discovery rate of roughly 
one new Local Group satellite per decade However, 
the Sloan Digital Sky Survey (SDSS) has been able to 
uncover a population of extremely low-luminosity satel- 
lite galaxies, which has roughly doubled the number of 
known satellites [3, 0, 0|- Determining how these new 
satellites fit in a given model for dark matter and cos- 
mology presents a very exciting theoretical challenge. 

The Cold Dark Matter (CDM) model predicts the ex- 
istence of hundreds of MW satellites that are expected to 
host galaxies at the faint end of the luminosity function 
Q . The ability of gas to cool and form stars in these low 
mass dark matter halos depends on a number of com- 
plex physical processes, such as supernova feedback, the 
photoionizing background, as well as mass loss due to 
tidal interactions Q. Despite the broad range of ob- 
served luminosities, the dark matter masses for all of the 
pre-SDSS satellites are constrained to within relatively 
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narrow range, approximately ~ [1 — 6] x 10 7 M Q within 
their inner 600 pc @, H| ■ Understanding this strong lumi- 
nosity bias at the low mass end is crucial to deciphering 
the formation of these dwarf spheroidal (dSph) galaxies 
as well as to constraining the nature of dark matter. 

In this paper we show that three new and nearby mem- 
bers of the Local Group discovered by the SDSS (Will- 
man 1, Coma Berenices, and Ursa Major II) are likely 
to have masses comparable to their more luminous coun- 
terparts. Initial estimates have already shown that these 
galaxies have mass-to-light ratios similar to or larger than 
the pre-SDSS dwarfs [iillQj]. With luminosities more than 
two orders of magnitude less than the pre-SDSS dwarfs, 
these new satellites are not only interesting in the con- 
text of galaxy formation at the lowest-mass scales, but 
also for indirect dark matter detection. The new dwarfs 
are very faint, but they contain large amounts of dark 
matter and are located quite nearby, making them ideal 
sites to search for signals of dark matter annihilation. 

Current and future observatories, including space- 
based experiments, such as GLAST [1JJ , as well as a suite 
of g round-based Cerenkov detectors, such as STACEE 
[II, HESS [3, MAGIC [3, VERITAS [H, CANGA- 
ROO 0, & HAWK [13, will search for the signal of 
7-rays from dark matter annihilations. The prospects 
for 7-ray detection from dark matter in well-known MW 
satellites with these observatories has been the subject 
of many previous studies [H, 0, H(| ■ All of these sys- 
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terns are interesting targets not only because of their 
large mass-to-light ratios, but also because they are ex- 
pected to have very low intrinsic 7-ray emission. This 
is in contrast to the situation at the Galactic center, 
where astrophysical backgrounds hinder the prospects of 
extracting the signal from dark matter annihilation [2lj |. 
Moreover, the known location of MW satellites makes 
a search of dark matter annihilation well-defined, unlike 
the search of completely dark substructure, which would 
rely on serendipitous discovery [H, [H, [H| • 

From the mass modeling of the dark matter halos, 
we provide the first determination of the 7-ray signal 
from dark matter from Ursa Major II, Willman 1, and 
Coma Berenices (Coma hereafter). These galaxies pro- 
vide promising targets for 7-ray detection for three rea- 
sons: 1) they are the among the closest dark matter dom- 
inated systems, 2) they are expected to be free from in- 
trinsic 7-ray emission, and 3) present data on their stel- 
lar kinematics suggest that their dark matter halos are 
as massive as the more well-known population of MW 
satellites. 

This paper is organized as follows. In section [TTJ, we re- 
view the theoretical modeling of dwarf dark matter halos 
and the calculation of the 7-ray flux. In section IIIII we 
present the likelihood function for determining the flux, 
and outline the theoretical priors in the modeling. In 
sections HVl and fVl we present the results and discussions. 



II. THEORETICAL MODELING 

When modeling the stellar distribution of a dwarf 
galaxy, it is important to determine the effect of ex- 
ternal tidal forces on the dynamics of the system. For 
the MW satellites we study, we can obtain an estimate 
of the external tidal force by comparing its magnitude 
to the internal gravitational force. The internal gravita- 
tional force is ~ a 2 /R and the external tidal force from 
the MW potential is ~ (220km/s) 2 i?/D 2 , where D is 
the distance to the dwarf from the center of the MW, 
and 220km/s is the MW rotation speed at the distance 
of the dwarf. The MW satellites are characterized by 
scale radii of R ~ 10 — 100 pc and velocity dispersions 
a ~ 5 — 10 km s _1 . For a typical distance of D ~ 40 kpc, 
the internal gravitational forces are thus larger by ~ 100. 

An additional estimate of tidal effects can be obtained 
by comparing the internal crossing times of the stars in 
the galaxy to the orbital time scale of the system in the 
external potential of the host. For a galaxy with scale 
radii and velocity dispersions given above, an estimate of 
the crossing time is given by R/a ~ 1 — 20 Myr. Assum- 
ing a rotation speed of ~ 220 km s _1 at the distance of 
these dSphs (~ 40 kpc), their orbital time scale in the 
MW potential is ~ Gyr. 

From the above estimates we conclude that it is highly 
unlikely that these galaxies are presently undergoing sig- 
nificant tidal stripping. These galaxies may have been 
tidally stripped before (for example, if their orbit took 



them closer to the center of the MW) , however the stel- 
lar core that has survived is faithfully tracing the local 
potential [H, [26|]. Of the galaxies we consider, only Ursa 
Major II shows strong evidence of past tidal interaction, 
as it is located on the same great circle as the Orphan 
Stream discovered by SDSS |3j, |4j . Thus we can proceed 
ahead with confidence in modeling the surviving stellar 
cores as systems in dynamical equilibrium. 

Line-of-sight velocities are widely used to determine 
the properties of the dark matter halos of dSphs (27L 
I28L [2^ | assuming spherical symmetry. For a system in 
dynamical equilibrium, the spherically symmetric Jeans 
equation gives the stellar line-of-sight velocity dispersion 
at a projected radius, R, from the center of the galaxy, 
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where the three-dimensional velocity dispersion, a r is ob- 
tained from 
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Here f3 is the stellar velocity anisotropy, and p+{r) is the 
density profile for the stellar distribution, which is ob- 
tained from the projected stellar distribution, I(R). The 
stellar distributions of the dSphs are typically fit with 
either Plummer or a King profiles [301 ] . The primary dif- 
ference between these fits is that Plummer profiles are de- 
scribed by a single parameter, rp, and fall-off as a power 
law in the outer regions of the galaxy, while King profiles 
are described by a core radius, rjc, and a cut-off radius, 
r cut , and fall-off exponentially in the outer regions. For 
the stellar distributions of Ursa Major II and Coma, we 
use the Plummer fits compiled in [l0(, and for Willman 
1 we use the King profile fit from [9(. These quantities, 
along with the distance to each galaxy and their respec- 
tive luminosities, are given in Table |TJ 

We model the distribution of dark matter in the dSphs 
with radial density profiles of the form 



p{r) 



Ps 



fT(l + f) 



S\3-7> 



(3) 



where p s is the characteristic density, f = r/r s , and r s 
is the scale radius. Numerical simulations bound 7 in 
the range [0.7 — 1.2], and the outer slope is ~ —3 [3l| . 
In order to compare the mass distribution in dSphs to 
dark matter halos in N-body simulations, it is useful to 



work in terms of the two parameters V max and 



the 



maximum circular velocity and the radius at which it is 
obtained. For example, for 7 = (0.8,1,1.2), r max /r s = 
(2.61,2.16, 1.72) . Thus, for a particular value of 7, the 
density profiles of dark matter halos can be described by 
either p s -r s or similarly by U max -r max . 

With the parameters of the halo density profile spec- 
ified, the 7-ray flux from dark matter annihilations is 
given by 
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dSph 


Distance (kpc) 


Luminosity (10 3 L©) 


Core Radius (kpc) 


Cut-off Radius (kpc) 


Number of stars 


Ursa Major II 


32 


2.8 


0.127 (P) 




20 


Coma Berenices 


44 


2.6 


0.064 (P) 




59 


Willman 1 


38 


0.9 


0.02 (K) 


0.08 (K) 


47 


Ursa Minor 


66 


290 


0.30 (K) 


1.50 (K) 


187 



TABLE I: The distance, luminosity, core and cut-off radii (Plummer (P), King (K)) for each of the dSphs we study. The last 
column gives the total number of stars used in the analysis. 



where rj± = D cos £ ± y r\ — D 2 sin 2 £, D is the distance 
to the galaxy, £ is the angular distance from the center of 
the galaxy, and r t is the tidal radius for the dark matter 
halo. Note that in the limit of D ^> r s , the flux scales 
as J p(r) 2 dr/D 2 , and in the particular case where 7=1, 
$ ~ p 2 r^/D 2 , with ~ 90% of the flux originating within 
r s - 

In Eq. ([4|), the properties of the dark matter particle 
are determined by 
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Here, E t h is a threshold energy, M x is the mass of the 
dark matter particle, (av) is the annihilation cross sec- 
tion, and the spectrum of the emitted 7-rays is given 
by dN/dE. Unless otherwise noted, we assume P rs 
10 _28 cm 3 s _2 GeV -2 , which corresponds to the most op- 
timistic supersymmetric dark matter models. However, 
we stress that the derived results can be rescaled to any 
dark matter model, by a simple rescaling of P. 



III. LIKELIHOOD FUNCTION AND PRIORS 

Observed stellar line-of-sight velocities place strong 
constraints on several important parameters describing 
the dark matter halos of dSphs. Two examples of these 
parameters are the halo mass and density at a character- 
istic halo radius, corresponding to about twice the King 
core radius (or about 600 pc for a typical dSph 0, [1]). 
More relevant to 7-rays is the quantity p 2 r 3 (the 7-ray 
luminosity is L ~ p 2 ^), which is typically determined 
to within a factor 3 — 6 with the line-of-sight veloci- 
ties of several hundred stars [ljj |2(j. The constraints 
on these parameters can be strengthened by including 
relations between similar parameters observed in numer- 
ical simulations. When combined with the observational 
constraints, these empirical relations in numerical simu- 
lations constitute a theoretical prior, delineating a pre- 
ferred region of the parameter space of the dark matter 
distribution in dSphs [20|, l32| . In this section, we discuss 
the implementation of this prior, and derive the general 
form of the likelihood function we use to constrain the 
7- flux from dark matter annihilations. 

We assume the line-of-sight velocities are drawn from 
a Gaussian distribution, centered on the true value of 



the mean velocity, u. This has been shown to be a good 
description of the well-studied dwarfs with line-of-sight 
velocities of several hundred stars [33[. Given the set of 
theoretical parameters, the probability to obtain the set 
of observed line-of-sight velocities, x, is 
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Here 9 is the set of parameters that describe the dSph, 
and the sum is over the observed total number of stars. 
The dispersion in the velocities thus has two sources: 1) 
the intrinsic dispersion, a t ,i(&), which is a function of the 
position of the i th star, and 2) the uncertainty stemming 
from the measurement, a m ^. 

We can simplify Eq. ^ by assuming that the mea- 
surement uncertainties are small relative to the intrinsic 
dispersion. This is a good approximation for well-studied 
dwarfs, which have intrinsic dispersions ~ lOkms^ 1 , and 
measurement errors ~ lkms -1 [HI]. Under this approx- 
imation Eq. ((6]) becomes 
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where the sum is now over the number of bins, Nb, for 
which the velocity dispersion is determined. The velocity 
dispersion in the i th bin is of j; and the number of stars 
in the i th bin is N % . We can use Eq. ([7]) if the observa- 
tions are given by line-of-sight velocity dispersions, and 
if the measurement errors are small in comparison to the 
intrinsic dispersion. This is the case for Ursa Minor, as 
discussed below. 

We describe the dark matter halos in terms of the pa- 
rameters 6 — (V max , r max , [3). We assume that [3 is con- 
stant as a function of radius, and let it vary over the 
range [—5 : 1]. We integrate Eqs. ^ and J7]) over these 
parameters and define the likelihood function for a fixed 
7-ray flux, /, as 

£(/) ex y"p(x|U max ,r max ,/3)F(K iax , ' max J 

x 8($(V max ,r ma , x ) - f)dV max dr max d/3. (8) 

Here we have assumed a uniform prior on f3, and the 
prior probability distribution for V^ nax ,r max is given by 
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dSph 


Mass < 100 pc [10 6 M ] 


Vmax [km S *] 


Ursa Major II 


o 1 +5.6 


qo+69 
zo -10 


Coma Berenices 


1 9+ 21 


19±f 


Willman 1 


1 3+ 15 


27-i 5 


Ursa Minor 


o q+1.9 
z.o_ 1 2 


30±i| 



10 



TABLE II: The masses within 100 pc and the maximum cir- 
cular velocities of the Milky Way satellites we study. Error 
bars indicate the 90% c.l. regions. No upper limit could be 
obtained for the maximum circular velocity of Willman 1. 



f (Knax, fmax)- This prior distribution is determined by 
the V ma x-?"max relation from CDM simulations. 

In order to determine P(V max , r max ), we need both its 
mean relation and its halo-to- halo scatter. For dark sub- 
halos that have been strongly affected by tidal interac- 
tions, the Vmax — i"max relation is strongly dependent on 
the nature of the potential of the host system, as these 
systems undergo varying amounts of mass loss as they 
evolve within the host halo. For example, Bullock and 
Johnston [35[ have embedded a disk potential in a MW 
size host halo, and found that the V max — r max relation 
of subhalos takes the form log(r max ) = 1.35[log(V max ) — 
1] — 0.196. We obtain a similar slope by examining the 
subhalos in the dark matter-only Via Lactea simulation 
of MW substructure [24| , however differences in the as- 
sumed cosmological parameters, and the absence of a disk 
potential in Via Lactea, translate into differences in the 



normalization of the V m 



relation. For Via Lactea 



subhalos, we find the normalization is reduced by ~ 30%, 
which implies reduced halo concentrations (larger r max 
for fixed Vnax)- 

We model the scatter in Vmax as a log-normal distribu- 
tion, with criogv.nax — 0.20. This provides a conservative 
estimate for the scatter in Vm ax as a function of r max 
for nearly the entire range of the subhalo mass function. 
At the extremely high end of the subhalo mass function 
(Vnax ^ 20kms~ 1 ), the scatter increases because in this 
range it is dominated by a small number of very massive 
systems that have been accreted into the host halo very 
recently. This increase in the scatter is simply because 
Via Lactea is only one realization of a substructure pop- 
ulation in a MW halo. We find that by excluding the 
extreme outliers in the Via Lactea mass function, the 
scatter is similar to the low mass regime. This is a sim- 
ilar result to those obtained in semi-analytic models of 
many realizations of the subhalo population (3ol [37j . 



IV. RESULTS 

A. Flux estimates for smooth dark matter 
distributions 

We now quantify the prospects for detecting 7-rays 
from dark matter annihilation in the three new dSphs. 
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FIG. 1: The 90% confidence level region in the p a — r a (top) 
Ps r t — ps (bottom) parameter space for Coma , Ursa Major II, 
Willman 1, and Ursa Minor. We marginalize over the velocity 
anisotropy and have assumed an inner slope of 7 = 1. The 
best-fit values are indicated with points. 



We first assume that the dark matter is distributed 
smoothly, and we discuss the implications of a boost fac- 
tor due to substructure in the next subsection. We use 
observations of these galaxies from the following refer- 
ences: Ursa Major II and Coma Willman 1 [38|. 
For these galaxies we have individual stellar velocities, 
so we use the likelihood function in Eq. ([5]). To make a 
connection to previous studies of dSphs, we compare the 
fluxes for these new dSphs to the flux from Ursa Minor, 
which is at a distance D — 66 kpc, and has a luminosity 
of L = 2.9 x 10 5 L©. Ursa Minor has the largest flux of 
any of the well-known dwarfs [20( . We describe the stel- 
lar distribution of Ursa Minor with a King profile, with 
tk = 0.30 kpc and r cut = 1.50 kpc (39|. For Ursa Minor 
we use the measured velocity dispersion from a sample of 
187 stars distributed evenly in 11 bins [39j, and we use 
the likelihood function in Eq. ([7]). 

In the top panel of Fig. |T]) , we show the 90% c.l. region 
in the p s — r s plane for each galaxy, where the best-fit 
values are denoted by points. Here we use the Vmax — '"max 
prior from section [TTT1 and we take the inner slope of the 
dark matter halo profile to be 7 = 1. In the bottom 
panel of Fig. {1} we show the 90% confidence level region 
p s plane. As seen in Fig. ([T]), the range of 
can take in each dSph is reduced with 
the inclusion of more stars (e.g. Ursa Minor vs any one 
of the other three dSphs). 

The constrained regions in Fig. ([T]) can be used to de- 
termine the masses, maximum circular velocities, and 7- 
ray flux probability distributions for each galaxy. In Ta- 
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FIG. 2: The probability distributions of the angular size sub- 
tended by r s for each galaxy. We marginalize over the velocity 
anisotropy and p s . The inner slope is fixed to 7 = 1. 



ble HH we show the masses within 100 pc and the max- 
imum circular velocities for each galaxy. The error bars 
indicate the 90% c.l. regions. 

To determine the flux distributions, we must first spec- 
ify a solid angle for integration. For optimal detection 
scenarios, the solid angle should encompass the region 
with the largest signal-to-noise. For the present work, 
we will integrate over a region where 90% of the flux 
originates. As discussed above, for the particular case 
where 7 = 1, 90% of the flux originates within r s . There- 
fore, in order to estimate the solid angle of integration, 
we have to first determine the maximum likelihood val- 
ues of r s . This is done by marginalizing over p s and (3 
with the Vmax — r max prior. The distributions of angular 
sizes are then obtained from /1 = tan -1 [r s /D], where D 
is the distance to the dSph. As shown in Fig. @, we 
find that given their similar size and roughly similar dis- 
tances, all three dSphs will emit 90% of their 7— ray flux 
within a region of ~ 0.2 degrees, centered on each dSph 
(for 7=1). Ursa Minor is the most physically extended 
galaxy, subtending the largest projected area on the sky. 

It is important to determine whether each of the galax- 
ies will be detected as point sources, or whether they will 
be resolved as extended objects. To determine this we 
compare their angular size to the angular resolution of 
7-ray telescopes. GLAST will have a single photon an- 
gular resolution of ~ 10 arcminutes for energies greater 
than 1 GeV, similar to the angular resolution of ground- 
based detectors (such as VERITAS) for energies greater 
than few tens of GeV. In the case where the detected 
number of photons is iV 7 > 1, the angular resolution of 
a detector is improved by a factor of 1/ y/N^. Therefore, 



FIG. 3: The probability distributions for the 7-ray fluxes from 
Coma, Ursa Major II, Willman 1, and Ursa Minor, marginal- 
izing over the velocity anisotropy, p a , and r a . We assume 
P = 10 -28 cm 3 s _1 GeV -2 and an inner slope of 7 = 1.0. 
We have assumed no boost from halo substructure, which in- 
creases these fluxes by a factor ~ 10 — 100. 



these galaxies can be resolved as extended objects, which 
in principle would allow a measured flux to determine the 
distribution of dark matter in the halo itself. 

Fig. ([3]) depicts the resulting flux probability distribu- 
tion for the three new dSphs and Ursa Minor. These 
are obtained by marginalizing over /?, p s , and r s and in- 
cluding the U max -r max prior. We set the inner slope to 
7 = 1, and integrate the flux over the solid angle that 
corresponds to 0.2 degrees from the center of the galaxy. 
We assume a value of P = P = 10~ 28 cm 3 s _1 GeV" 2 , 
but the result can be scaled to any dark matter candi- 
date with a different value of P by simply multiplying 
the flux distribution by a factor of P/Pq- 

The relative proximity of the three new dSphs, and 
their comparable sizes, results in 7-ray fluxes that are 
roughly similar. For P ss Pq, the likelihood peaks at ap- 
proximately $0 ~ 10 cm s , with a spread of nearly 
an order of magnitude. Thus Ursa Major II, Coma, and 
Ursa Minor all have comparable fluxes, and Willman 1 
has a most likely flux that is about three times larger 
than Ursa Major II or Coma. 



B. The effects of the inner slope and substructure 
boost factors 

Understanding the distribution of dark matter in the 
inner regions of the dSphs also has important implica- 
tions for detection of a 7-ray flux. However, when vary- 
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Flux [photons cm~ a s~'] Flux [photons cm" 8 s~'] 

FIG. 4: The 7-ray flux probability distributions for Coma and 
Ursa Minor for inner slopes of 0.8 (long-dashed), 1.0 (solid), 
and 1.2 (short-dashed). We marginalize over the same quan- 
tities as in Fig. ([3j . The value of P is the same as in Fig. ([3J . 



ing the inner slope, we must also be certain to vary all 
of the other halo parameters so as to remain consistent 
with the line-of-sight observations. In order to quantify 
the effects of varying the inner slope, we marginalize over 

Kiax 5 ''max 

and (3 for profiles with different values of 7. 
In Fig. ([4]) we show the effects of varying 7 for Coma 
and Ursa Minor. The shifts in the flux distribution func- 
tion are not only a result of varying the inner slope but 
also come from the constraints imposed by the data on 
the density profile parameters p s and r s . The relative 
amount of the shifts can be understood by considering 
the best-fitting values of r s (where the majority of the 
7-ray flux comes from) as compared to the core radii of 
the systems. When the core radius is similar to r s , the 
shifts are larger for varying 7, as in the case of Ursa Mi- 
nor. However, when the core radius is much smaller than 
the fitted values of r s , variations in the inner slope are 
less significant than the change induced by p s , as is the 
case for Coma. 

The presence of substructure in dark matter halos is 
firmly established on theoretical and numerical grounds. 
Dark matter halos are approximately self-similar, and 
substructure is expected to be present in all dark mat- 
ter halos with mass greater than the cut-off scale in the 
primordial power spectrum, set by the kinetic decoupling 
temperature of the dark matter particle (for detecting the 
smallest dark matter halos in the Milky Way, see [23|]). It 
is therefore natural to expect that these galaxies contain 
substructure if they consist of CDM. 

The density enhancement over the smooth distribution 
of dark matter leads to an enhancement in the total an- 
nihilation rate, typically quantified in terms of a "boost" 
factor. As was shown in [2fJ, the boost factor cannot at- 
tain arbitrarily large values, but instead is bounded to be 
less than ~ 100, with the exact value depending on the 
cut-off scale in the CDM mass function. The boost is a 
multiplicative quantity, so the effect of dark substructure 
is simply accounted for by scaling the fluxes in Fig. ([3]) 



by the appropriate boost factor. 



C. Detection prospects 

As is shown in Fig. ([3]), the flux probability distribu- 
tion functions peak around $0 ~ 10 _10 cm _2 s _1 with- 
out including any enhancement to the signal from sub- 
structure. Here, we will assume a conservative value 
of 10 for the boost factor, and discuss the prospects 
for detecting the three new dSphs with 7-ray instru- 
ments. We can make simple estimates for the likeli- 
hood for detection by adopting the specifications of par- 
ticular 7-ray detectors. We will use two examples: a 
space-based experiment, GLAST, and a ground-based 
Cerenkov telescope, VERITAS. For GLAST, if we as- 
sume an orbit-averaged effective area of A c g w 2 x 10 3 cm 2 
and an exposure time of t exp = 10 years, their product 
is Bq = A c fft cxp 3 x 10 11 cm 2 s. A 50 hour exposure 
with VERITAS (A eS » 10 8 cm 2 ) has B v «2x 10 13 cm 2 s. 
Naively, for a fixed value of P, a ground-based detector 
seems more sensitive because By > Bq. However, the 
backgrounds for a ground-based detector are also larger 
and include a component from the hadronization of cos- 
mic rays in the atmosphere. 

As an example, for a fixed value of P = Pq, and a 
solid angle that corresponds to an angular size of 0.2 
degrees, the number of photons detected by GLAST is 
Nj,g ~ 300. The dominant source of background for 
GLAST is the Galactic diffuse emission (dNn/dE = 
1.2xlO- 6 [Sth/GeV]- 11 cm- 2 s- 1 sr- 1 GeV" 1 40]). If we 
assume an energy threshold of E t h ~ 1 GeV, then the 
number of background photons is Nb — Bq cL/Vb / dE « 
250, which means that the new satellites will be de- 
tected at approximately a N^^q/WNb + N 7i g ~ 12c 
level. A similar estimate can be obtained for VERI- 
TAS. The number of photons detected above 50 GeV 
in an instrument with an effective area times exposure 
~ By is N 7 w 2 x 10 4 . The dominant contribution 
to the background are photons that originate from neu- 
tral pion decays from the nuclear interactions of cosmic 
rays in the upper layers of the atmosphere (dNs /dE = 
3.8xl0- 3 [ J B t h/GeV]- 2 - 75 cm- 2 s^ 1 sr- 1 GeV" 1 'Al}). The 
number of background photons from pion decays is ap- 
proximately Nb ~ 2 x 10 7 , therefore the 3 dSphs could 
be detected at a v « 5er level. Understanding and 
discriminating against the background contamination in 
Cerenkov telescopes is very important in improving the 
prospects for detecting dSphs of the Milky Way. 

As shown in [20( , the large number of stellar velocities 
obtained in the older dSphs allow useful 7-ray flux ratios 
between different dSphs to be determined. For the three 
new dSphs considered in this work, the kinematic data is 
not good enough to play the same game. Clearly, more 
stellar velocities will shrink the allowed region of p 2 r 3 
parameter permitting robust estimates of flux ratios be- 
tween the galaxies studied here and the rest of the Milky 
Way satellites. 
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V. DISCUSSION & CONCLUSIONS 

In this paper, we have modeled the dark matter 
distribution in three recently-discovered Milky Way 
(MW) satellites (Ursa Major II, Willman 1, and Coma 
Berenices), and have presented the prospects for detect- 
ing 7-rays from dark matter annihilations in their halos. 
We show that the expected flux from these galaxies is 
larger than the flux from any of the higher luminosity, 
more well-known (pre-SDSS) dwarfs. There are two rea- 
sons for this surprising result: 1) the masses of these new 
dwarfs within their stellar distributions are similar to the 
masses of the well-known, larger luminosity dwarfs, and 
2) all three new galaxies are closer than the other well- 
known dwarfs. The implied mass-to-light ratios, ~ 1000, 
of these new dwarfs makes them the most dark matter 
dominated galaxies in the Universe. 

Our estimates show that it is unlikely that the ob- 
served stellar distributions are presently undergoing tidal 
disruption. However, this does not mean that they have 
been free from tidal interactions in the past, but rather 
that the surviving stellar core can be faithfully modeled 
as a system in dynamical equilibrium. By including the 
VJnax — ''max CDM prior, we have naturally accounted for 
tidal effects in the mass modeling, since this V max — ''max 
relation in fact comes from dark matter halos that have 
experienced tidal stripping. 

One of the galaxies we consider, Ursa Major II, may 
be a candidate for past tidal disruption, given that it 
is positioned on the same great circle as the Orphan 
Stream of stars, which was also recently detected by the 
SDSS @, i]. This is consistent with the findings of Si- 
mon and Geha flO] , who have recently investigated the 
possibility of tidal disruption in Ursa Major II, as well as 
all of the other new dwarfs, using proxies such as gradi- 
ents in the observed velocity distribution and metallicity 
of the stellar populations. Given the total mass-to- light 
ratio we have determined for Ursa Major II, tidal strip- 
ping will have only been significant if its pericenter is ~ 
3 times closer than its present distance. Future observa- 
tions of the stellar distributions in Ursa Major II, and all 
of the other new faint dwarfs, will be important in deter- 
mining bound and unbound stellar populations. With a 
larger sample of stars from a galaxy such as Ursa Major 
II, it will be possible to remove unbound and interloping 
stars with techniques similar to those presented in [251 ] - 
Upon removal of stars unbound to the galaxy, these au- 
thors show that in most cases the true bound mass of the 
system can be recovered to typically better than 25%. 

As a very conservative check for the effects of member- 
ship uncertainties, we have redone the analysis for each 
of the galaxies by just keeping the stars within the inner 
half of each galaxy, where the surface densities are the 
most well-determined. We find that the peak of the flux 
likelihood is shifted by a small amount relative to the 
l-er widths in Fig. ((3]). Note however, that even when 
including the entire population of stars in the observed 
samples, in all cases equilibrium models provide adequate 



descriptions of the dynamics of each system. 

Unresolved binary star systems also introduce a sys- 
tematic that may effect the fluxes we have presented. 
Olszewski, Pryor, and Armandroff [42| have determined 
the effect of binaries on the velocity dispersion of two 
of the most luminous dwarfs, Draco and Ursa Minor, by 
inferring the binary population of these systems using 
multiple epoch observations. They find a velocity dis- 
persion of ~ 1.5 km s^ 1 due to binaries, and a probabil- 
ity of 5% that binaries elevate the velocity dispersion to 
4kms~ 1 , which is still less than the velocity dispersion 
of the three new dwarfs. Thus if Ursa Major II, Willman 
1, and Coma have binary fractions similar to Draco and 
Ursa Minor, their observed velocity dispersions are not 
significantly affected by binaries. We note that this is 
consistent with recent estimates of the binary fraction in 
low density Galactic globular clusters [HI]. 

A strong test for the presence of binaries is to examine 
the distribution of measured velocities. The velocity dis- 
tribution due to internal motion in binary systems should 
be flat due to the observed broadness of the period distri- 
bution of binaries [3] . The observed period distributions 
depend on spectral type and age of the system (among 
other variables) but are all broad with dispersions of 
about two orders of magnitude, which seems to be con- 
sistent with theoretical expectations (45|. Further, there 
is no observational evidence or theoretical argument that 
suggests that the period distribution should be sharply 
suppressed for all periods below ~ 1000 years (roughly 
velocities larger than ~ 3 km/s). Thus if there is a large 
contribution from internal motion in binary stellar sys- 
tems to the intrinsic velocity dispersion of these dwarfs, 
we expect to see a significant tail of high velocities. This 
is not observed and hence we can be confident that mea- 
sured velocity dispersion is tracing the total mass in the 
dwarf galaxy. 

The new, ultra-low luminosity galaxies represent an 
interesting confluence of astronomical and 7-ray studies. 
Future kinematic studies of all of these new dwarfs will 
further reduce the errors on the mass distributions, and 
sharpen the predictions for 7-ray observatories searching 
for signatures of dark matter annihilations. 
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